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Abstract 

An extension of the left-right symmetric model has been con- 
structed which gives in a natural way the three lepton decay modes of 
the proton which have been suggested as an explanation for the atmo- 
spheric neutrino anomaly. We write down the potential which after 
minimization gives the proper choice of the Higgs spectrum. With 
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this Higgs spectrum we then study the evolution of the gauge cou- 
pling constants and point out that for consistency one has to include 
effects of gravity. 
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1 Introduction 



One expects to see produced in the atmosphere twice as many muon neutrinos 
as electron neutrinos since detectors cannot distinguish between neutrinos 
and antineutrinos. The two water-Cerenkov detectors give a result which is 
a factor of two smaller for the ratio R = Niy^/Niy^), a ratio in which many 
systematic uncertainties are expected to cancel. The results are 

Robs/Rue = 0.60 ± 0.07 ± 0.05 

from the Kamiokande experiment [p]] (based on 6.1 Kton year) and 

Robs/Rue = 0.54 ± 0.05 ± 0.12 

from 1MB |2j (based on 7.7 Kton Year). The experiments look for "contained" 
single prong events which are caused by neutrinos with energies below 2 GeV. 
The ratio R is estimated from the relative rates of sharply defined single rings 
(muon-like) and diffused single rings (electron-like). (Other reported values 
§ for this ratio are: Frejus [g 0.87 ± 0.21 (1.56 Kton Year), NUSEX § 
0.99 ± 0.40 (0.4 Kton year), SOUDAN II § 0.69 ± 0.19 (1 Kton year)). 

Although this atmospheric neutrino anomaly has a popular explanation 
within the neutrino oscillation framework 0, there is an alternative expla- 
nation based on three lepton decays of the proton || . The single ring events 
have been analysed within the proton decay interpretation where it is argued 
that, if the proton decays into a positron and two neutrinos with a lifetime 
of r(P — > e + vv) ~4x 10 31 years, then the excess observed electron events 
could be due to proton decay events |§. The lifetime for this particular decay 
mode of the proton || is consistent with the present limit |10[ for the ex- 



pected dominant decay mode of the proton r(P — > e + 7r°) > 5 x 10 32 yr. The 
possibility that this particular decay mode might dominate over other usual 



decay modes was considered earlier on general grounds |TT], O, [13| where it 
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was pointed out that it is difficult to have light neutrinos in the final decay 
product. 

In most models the left-handed neutrinos are light and the right handed- 
neutrinos are heavy. Thus the decay modes are restricted to 

P — ► e + v L u L or P —> e + v L c v L or P — > e + v L c i' L c . (1) 



Recently [T3| it has been pointed out that these decay modes are allowed 
in the framework of certain left-right symmetric models. In this paper we 
construct an explicit model in which three lepton decay mode of the proton 
is the dominant one. We write down the general form of the potential and 
show that the minima of the potential are consistent with the choice of Higgs 
structure. We then study the evolution of the gauge coupling constants 
including non-renormalizable interactions which may arise from the Planck 
scale physics |14|, [15 . 
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2 A Proton Decay Mechanism 



We work in the framework of the left-right symmetric model [IT], [12, |l^] and 
start with the symmetry breaking chain 

SU(4) x SU{2) L x SU{2) R [= G PS ] 
SU(3) C x SU(2) L x SU(2) R x 17(1) (B -x) [= G LR ] 

M R 



SU(3) C x SU(2) L x C/(l) y [= G 



std\ 



Mi 



Stf(3) c x [/(l) f 



In the minimal left -right symmetric model [II, 12, |16[ the Higgs scalars 
consist of the following fields. The group Gps is broken by the vacuum 
expectation value (vev) of the field H which transforms as (15,1,1) under 
the group Gps- The right handed group is broken by the vev of a right 
handed triplet Higgs field A R = (1,1,3,-2) C (10,1,3). By left-right parity 
this will imply the existence of the left handed triplet field = (1,3,1, 
-2) C (10,3,1), which gives Majorana mass to the left-handed neutrinos and 
whose vev should be < 1 GeV. (Where there are four numbers, the first three 
correspond to the representations of SU(3) C x SU{2) L x SU{2) R while the 
last shows the U(l) quantum numbers). 

Finally the electroweak symmetry breaking takes place through the vev 
of a doublet scalar field = (1,2,2,0) C (1,2,2). This field (j> a ^ so gives 
masses to the fermions. However this does not reproduce the right quark- 
lepton mass ratios. For the right magnitude of the quark-lepton mass ratios 
we require yet another field f = (1,2,2,0) C (15,2,2) 0. The SU(3) C 
singlet component of this field £, which acquires vev, has different Clebsch- 
Gordon coefficients for the SU(3) C and the U(l) part of S'f/(4). Hence they 
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contribute to the quark and lepton masses with different coefficients. As a 
result, suitable combinations of and £ can reproduce the right quark-lepton 
mass ratios. 

It was pointed out in ref |H| that with this minimal scalar content it is 
possible to get the decay mode required to explain the atmospheric neutrino 
problem. For this we need the SU(3) C color triplet components of the fields 
Al and £, which we represent by A| and £ 3 respectively. (Al is (10,3,1) 
under Gps and the (1,3,1,-2) component acquires a vev; the 10 represen- 
tation of S77(4) decomposes under SU(3) C as 6 + 3 +1. Similarly, the 15 
representation of S77(4) decomposes under SU(3) C as 8 + 3 + 3 + 1). Then 
the Yukawa couplings, 



C Y uk = f q i(^l L )&T + fdi{d R c l L c )e (2) 
and the quartic scalar coupling, 

C s = A^AlAl^ 1 (3) 

give the (B—L) conserving proton decay P — > e^viy^ through the diagram 
of figure 1. This diagram will also give, with equal probability, the decay 
mode, P — > /ii + z/ L z/ L c . 

Such a proton decay mechanism will give equal number of sharp single 
rings (muon-like events) and diffuse single rings (electron-like events). Since 
the proton decay events contribute to both electron- and muon-like events 
this seems to imply that the reduction of the ratio R cannot be explained 
by proton decay events. However, the weighted average of the two processes 
with ratios R(proton decay) = 1 and R(atmospheric neutrino) = 2 (the 
theoretical expected ratio for the muon-to-electron events if they have their 
origin only from the atmospheric neutrinos) can in fact explain the atmo- 
spheric neutrino anomaly. 

To see this we note the observed numbers of electron-like [muon-like] 
events n e (obs) [n^obs]] are the sum of the numbers of electrons [muons] pro- 



6 



duced by the atmospheric electron [muon] neutrinos v e through scattering 
inside the detector n e (atm) [n^atm)] and from the decays of the protons into 
g + VlVl c \P' Jr VL^L c \ inside the detectors n e (prot) [n^prot)}. That is, 
n^obs) n^atm) + n^prot) 

Robs = — 7-tt = — rz — r~; — 7 — iv ~ u.b/tMc ~ ±2. 

n e {oos) n e {atm) + n e {prot) 

In ref f| it was assumed that the proton decays into e^Ujjujf (and not 
muons), i.e., n^{prot) = 0. They found that, by doing a Monte Carlo simu- 
lation to obtain the proton lifetime, the atmospheric neutrino anomaly could 
be achieved with a proton lifetime of t p ~ 4 x 10 31 years. In the above 
relation this corresponds to n e (prot) ~ (2/3)n e (atm) ~ (1/3)71^ (aim). 

In the present scenario the proton decays into both electrons and muons 
so that n e (prot) ~ n^(prot). Thus for the explanation of the atmospheric 
neutrino anomaly we require n e (prot) ~ 4n e (atm) ~ 2n^{atm). Since the 
number of proton decays is increased to give the same R b s there will be 
a reduction in the proton lifetime by a factor of 6. Thus in this scenario 
we can explain the atmospheric neutrino anomaly with a proton lifetime 
t p ~ (2/3) x 10 31 years, which is still consistent with present experiments on 
proton decay. 

The amplitude for the process is given by 

A = XPTf f fd f\ (4) 

where, (t; 1 ) = ((f)) = 250 GeV, X pr is the strength of the quartic coupling 
defined in Eq. (|3|) and the f q i, fdi are the Yukawa coupling constants. 

Then, taking reasonable values for the quartic and the quadratic Yukawa 
coupling parameters, X pr ~ and / ~ 10~ 3 , say, requires m^3 and m&3 
to be relatively light. For the proton decay mode P — > e i + z/ i z/ i c to be 
10 31 years to explain the atmospheric neutrino anomaly, it has been argued 



IT], [L2| that the mass can be as light as about a TeV, which requires 



m&3 ~ few TeV. This can also be achieved naturally |L3 
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3 A General Left— Right Symmetric Potential 



We now concentrate on the masses m^, and . In earlier references [JTTJ, [12 
13fl two new mechanisms were proposed which could give rise to appropriate 
masses for these fields. Here we check the consistency of these two different 
mechanisms when the complete potential with all the scalar fields is written 
and minimized. First we describe the two mechanisms which keep these two 
fields m^3 and light. 

For light mg3 it was argued that if there exists a field £' = (15, 2, 2) 



which can mix with the field £, then fine tuning can give a large mass to one 
combination of the fields £ and £' and keep the other combination with a 
light mass. However, this will also keep the masses of the color octet and the 
color singlet light, which is undesirable from the point of view of evolution 
of the gauge coupling constants. This problem is avoided if instead of £' 



we introduce a field x = (6, 2, 2) under Gp,s- If the symmetry group Gps is 
embedded in the unified group SO (10) then this field is contained in a 54 
plet of 5*0(10), which is required to break the symmetry of the large group. 
The mixing of the field \ with £ can then give a mass matrix which may be 
fine tuned to give only a light color triplet field. We shall discuss the details 
of this mechanism at a later stage. In the rest of the article we shall use the 
field x an d not 

For the field A 3 to remain light we have to alter the way in which the 
left-right symmetric model gets broken, although, as we shall see later, this 
particular method will make the mechanism in the last paragraph consistent 
with the minimization of the general potential. For this purpose we introduce 
a singlet field r\ = (1,1,1,0) C (1,1,1), to break the left-right parity (usually 
this is referred to as D-parity) at a different scale from the left-right symme- 
try breaking scale Mr WJ§. This field transforms under D as rj — > —rj. The 
scalar and the fermionic fields transform under D-parity as Al^r — » Ar^l 
and if>L,R ~^ ipR,Li while and £ stay the same. Then with the field rj we can 
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add new terms to the lagrangian, 

C v a = -M vV (A R A R - A{A L ) - X v r] 2 (A[A L + A R A R ) (5) 

In theories where the triplet Higgs breaks D-parity along with SU(2) R 
we have m AL = m AR . The masses of the fields A L and A R are not be the 
same when D-parity is broken by the vev of the field 77. When 77 gets a 
non-zero vacuum expectation value, the masses are given by, 

m A L = m A~ M v(v) + and m\ R = m\ + M^rj) + \ v (r]) 2 

where wa is the mass at which the left -right symmetry breaking of SU(2) R 
is broken spontaneously. So, in the absence of the r\ field, both these A fields 
will have mass ~ m A . With < 77 > present the parameters in the three terms 
can be tuned to make m& L vanish. The field will then acquire mass of 
the order of a TeV from radiative corrections. The same sets of parameter 
will also make m AR heavy and lead to a solution 

(77) ~ (A R ) > (A L ) and M v « m AR « (A R ) « m A and tji Al < (A R ) (6) 

Thus we can have m A 3 ~ m A L ~ f ew TeV even when m AR ~ (A R ) is as 
large as 10 10 GeV. 

We now write the most general potential with all the fields present in the 
minimal left-right symmetric model with the additional fields £, x an d the 
D-parity odd-singlet field 77. We then show the two mechanisms required 
to keep the color triplet fields light are consistent with the minima of the 
potential. To simplify the expression we define, 

0i = ; 02 = r 2 (f)* 1 T2 ; £1 = £ ; £ 2 = t 2 £ 1 t 2 

The most general potential with all the fields is, 

\/(0 1 ,0 2 ,A L ,A i? ,6,6,^x) = V^ + V A + V v + Vs + V x 

+v v4> + V vA + V A4> + V# + V A £ + V v z + V X £ 

(7) 
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where the different terms in this expression are given by, 

i,j i,j,k,l 

+ E K 3 ki trtfUAh) 
V A = V MA{A L ) + tr(A R A R )} + Pl [tr{A{A L f + tr(A R A R ) 2 } 

+p 2 [tr(AiA L A{A L )+tr(A^ R A R A^A R )]+p 3 tr(A{A L A R A R ) 

V x = M 2 x tr( X ] X) + A^tr(x f x)f + AMx f XX 1 x) 

Vk = +E«u [tr(AlA L ) +tr(Al ? A i? )] fr(0i&) + Eft* [ ^(A^A^^J) 
+tr(A^A fl 0}^.)] MA^A^t) 

V vA = -Mr, V [tr(A[A L ) - tr(A R A R )} + \ r, 2 [tr(A^A L ) + tr(A R A R )] 

V v4> = Y, 6 ij V 2 tr{4(pj) 

i,j,k,l i,j,k,l 
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+ £ay [ tr(A{A L )+tr(A R A R )\ tr(^) 
+ [ tr(A{A L U}) + ir(AW£te)] 



^ = P r/[tr(exA fl ) - tr{tx& L )] + M [ir(&A fl ) + tr(&A L )] 



We have not written the 577(4) indices explicitly for simplicity. For exam- 
ple, if we include the 577(4) index, the term p 2 tr (A^A^A^A^) in our nota- 
tion will actually mean two terms, p%tr (A[ q A Lq A^A L/3 ) and p b 2 tr(A^ A Lp Af A La ). 
However, as far as the minimization and the consistency of the model is con- 
cerned, we only have to replace p 2 by (p% + p\). Otherwise the rest of the 
analysis will be unaltered. A more detailed analysis with explicit SU(A) 
indices will not constrain or relax any of the constraints in this model. 

The vacuum expectation values (vev) of the fields have the following form: 



< <p > = 

< 4>> = 

<£> = 

< i] > 



k 

k' 

k! 

k 

k 

k' 



< A L >-. 



< A R >= 



= Vo 



<Z>- 
< X >= o 




v L 

v R 
k' 
k 



The notation needs some clarification. For the fields <fi and £ we have 
used the representation in which rows correspond to the SU(2)l quantum 
numbers (+|, — |) and columns correspond to the SU (2)r quantum numbers 
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(— |, +2)- The field is a singlet under the group £{7(4) and hence it 
has a £{7(4) matrix representation diag(l, 1, 1, 1). On the other hand the 
field £ transforms under £{7(4) as a 15 representation. Under the SU(3) 
subgroup of £{7(4) the 15 decomposes as 8 + 3 + 3 + 1. The £{7(4) matrix 
representation of the singlet is a traceless diagonal matrix which is a unit 
matrix in the SU(3) space. Hence the SU(A) matrix representation of the 
component of £, which is a singlet under both the SU(3)c and the {7(1) 
subgroups of £{7(4) and which acquires a vev, is diag(l, 1, 1, —3). The SU(2) 
representations are as above. Thus these fields <fi and £ contribute differently 
to the quarks and leptons masses, and hence a proper combination of the 
two fields give the correct mass relations between the quarks and leptons 
[fn}| . For the fields At and A# we used the 2x2 triplet representations of 
SU(2). Thus the components r 1 — i r 2 , which has the isospin +1 and hence 
charge neutral, acquire vev. (The electric charge is T 3 l + T 3 r + (B — L)/2. 
For these representations B — L = —2. So the charge neutral component 
should have T 3 ^ or T 3 ^ = +1, meaning they should contract with r 1 —i r 2 ). 
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4 Minimization of the Potential 



It is almost impossible to minimize the potential with respect to all the fields 
and then find the absolute minima. For this purpose one needs to simplify 
the problem considerably. As a first approximation one can extremize the 
potential with respect to all the fields and then substituting the vevs of the 
different fields to check if there is any inconsistency. In the minimal left-right 
symmetric potential, i.e., without the field rj, there are no linear terms in 
any field so the usual practice is to replace the various fields by their vevs 
then extremize it with respect to these vevs and finally check for consistency. 
We shall also follow the same procedure, but we need to take care of the 
extra linear terms present in the potential. For these linear terms, we shall 
afterwards minimize the potential with respect to those fields which do not 
acquire any vev. The vanishing of these derivatives after substituting for the 
vevs will then impose new constraints which also have to be satisfied. 

After the spontaneous symmetry breaking, when the fields acquire a vev, 
the potential contains terms with k, k', k, k , vl and v R . We need only terms 
involving vl and v R . These are given by, 

V = -» 2 (vl + vl) + ^(vi + v* R ) + ^(vlv 2 R ) + 2v L v R [( ln 

+l22)kk' + 7l2 (/c 2 + k' 2 )} + {v 2 L + v 2 R ) [(a n + a 22 + (3 U ) k 2 
+(a u + a 22 + P22) k' 2 + (4ai2 + 2p 12 ) kk'\ 

-Mr, Vo {v\ - V R ) + \ T] 2 (V 2 L + V R ) 

+ 4) [( fl n + a 22 + b n + Ki) k 2 + (a n + a 22 + b 22 + \ p 22 ) k' 2 
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+ (4a 12 + b 12 + A? 2 r ) ~kk>] + 2v L v R [(c n + c 22 + A£ + A£)fc& 



+ (c 12 + A^)(P + fc' 2 )] (8) 

where we have defined the new parameters p = 4(pi + p 2 ) and p' = 2p%. 

The minimization of this potential gives a constraint on v^ and vr. In- 
stead of minimizing this potential with respect to the fields vl and vr sepa- 
rately, we consider a combination, J^wr — J^-Ul = 0, which gives a relation 
among the fields vl and vr 

Pik 2 + (3 2 k 2 fQ . 
VlVr = \, j (9) 

where, Pi = 1^\ 2 \ and (3 2 = 2(ci2 + Ai 2 ) and we assumed k' « k and k « k. 
This allows us to have a very tiny vev for the left-handed triplet field Ax 
while keeping the vev of the right-handed triplet field Ar very large. This 
is in agreement with what we required in Eq. (^|). 

Now consider the linear terms involving x given by V x ^. These terms allow 
for the correct mixing between the color triplet components of the fields £ 
and x- The mass matrix for £ 3 and x 3 is now given by 

M <\) <"» 

where, a = m 2 , d = M 2 and b = c = (Pr] + M)vr, and where we assumed 
all rriijS are equal to m. If we now fine tune parameters to make det.M = 0, 
i.e., 

(mM x ) 2 = (P Vo + Mfv 2 R 

then one of the mass eigenvalues is zero. This fine tuning requires that 
M must be negative and (Pt]o + M) be very small and negative. In fact, 
\Pr]o + M\ has to be of the order of ~ This massless field will get a 

mass of the order of ~ TeV after radiative corrections during the electroweak 
symmetry breaking are included. 
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The terms linear in x required for this mechanism have another effect 
which was not transparent when we did the minimization of the potential 
with respect to the vevs of the various fields. If we first minimize with 
respect to the field x an d then substitute for the vevs of the various fields 
(which is not usually done since that complicates the calculation), then there 
is an additional constraint, 

Pyo + M 



M 2 

In other words, for Vl ~ ~~m^-> we nee d 



This is consistent 

N1 R 



This means that to satisfy Eq. (§), we require \Pr]o + M \ << \Pi] — M\. 

Pyo+M 
Pvo-M 

with the fine tuning used to keep the color triplet component of £ light, for 

M 2 

which \Pr] + M\ ~ and |P?7o — M\ ~ Mr. This would not have been 
possible if the field 77 were not present. For example, in the original paper 
|T^j where the field x was introduced and the field 77 was not required, this 



method would have led to an inconsistency. In the absence of the field 77 
minimization of the potential with respect to the field x would have given 
a constraint v^ = —Vr, which is inconsistent with the LEP data. Here we 
need the field 77 to keep the left-handed color triplet light and the fine tuning 
makes it consistent. 

We now turn to the question of light A 3 . As mentioned earlier, to have 
Ma ~ few TeV, we require the coupling of rj and A as in Eq. (|^). Also, 
in the most general potential the only term which contributes at the level 
of fj ~ M& R is V v a, which is exactly the same as in (Eq. H). Thus the 
field A 3 ~ Al can be massless at that level. Then during the electroweak 
phase transition this will again acquire mass through radiative corrections 
of the order of a few TeV. This is more natural in supersymmetric theories 
where the radiative corrections induce mass of the order of supersymmetry 
breaking scale, which are usually of the order of a few TeV. 

The mechanism just mentioned to make A 3 light has one drawback. It 
makes the other components of A^ also very light. For example, the A 6 can 
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now mediate n — n oscillation, which has to be suppressed. This problem is 
similar to the doublet-triplet splitting of any other grand unified models. In 
the present scenario we assume that although this field is light the coupling 
of A 6 is very small, which can make the model safe. However, this is not 
the best choice and if one can find some good solution to the doublet-triplet 
splitting in other GUTs, then one has to incorporate the same mechanism 
here in future. With our present assumption that the Yukawa coupling of A 6 
is very small, we now have to check the consequences of these fields in the 
evolution of the gauge coupling constants. 

Finally, we point out that the quartic coupling (given by Va$) required 
for generating the required diagram is also present in the general potential. 
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5 Evolution of the Coupling Constants 



There have to be many light scalars for the present scenario to work. These 
scalars may destabilize the unification of the gauge coupling constant at 
the unification scale. In the evolution of the gauge coupling constant with 
these Higgs scalars included and with the mass scales as above it is im- 
possible to have unification of the gauge coupling constants using the LEP 



constraints |L8| on sin w and a s . However higher dimensional operators, which 
might originate from Planck scale physics such as quantum gravity or com- 
pactification of Kaluza-Klein theories or Superstring theories, can save the 
situation |L4], |15| . Thus if the three lepton decay mode of the proton survives 
all of the experimental tests, then we may have an indication that Planck 
scale physics is actually modifying the boundary conditions of the gauge 
coupling constants near the unification scale. 

In our analysis we include the effect of the non-renormalizable terms 
arising from Planck scale physics from the beginning using the notation and 
method of reference We write down the generalized renormalization 

group equations in which the Planck scale effects are parametrized in terms 
of four extra parameters. We recover the usual relations between the coupling 
constants in the absence of gravity by setting the extra parameters to zero. 
It is when we do this we obtain a contradiction and the equations fail to 
provide unification. 

The evolution of gauge coupling constants with the energy scale is gov- 
erned by Renormalization Group Equations (RGE). Here we consider the 
RGE in one loop approximation i.e. the gauge fields fermionic fields and the 
scalar fields contribute to the evolution of the gauge couplings via one loop 
graphs only [19|]. In this approximation the renormalization group equation 
takes the form, 

V d = 26iO:- 0) 
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where, = |^ and the beta function is given in the following generic form. 



1 

Arc 



11 »r 4 l m 

3 3 ; 6 



(12) 



Here N = 1,2, 3 or 4, the number of neutrinos is always 3 and the scalars 
take on the values discussed below. Since there are a large number of scalar 
fields present in our model their contribution will be substantial despite the 
suppression by a factor of 6 in the beta function of the scalar term (T s ). We 
list the scalar fields that contribute to the RGE at different energy scales in 
Table. [I]. For simplification of notation we write M c for Mps in this section. 

We embed the Pati-Salam group Gps into a larger GUT group SO (10). 
The 5*0(10) symmetry is broken by a 54-plet of Higgs field S at the scale 
Mxj. The E is a traceless symmetric field of the S'O(IO) and the vevs of S 
which mediates this symmetry breaking are given by, 

<E> = 7/^ S diag(l, 1, 1, 1, 1, 1, -| -| -| -\). (13) 



where, S = y -^^My and «g = 9$/^ is the GUT coupling. The vev of a 
45-plet field H breaks the symmetry group Gps to Glr, 



1 



O33 I33 O34 \ 
(H) = -^=-. H I -I33 33 34 (14) 

1 o 43 o 43 o 44 y 

where, mn is a m x n null matrix and l mm is a m x m unit matrix. 

The (1,1,1) component of the 54-plet field £ breaks the SO (10) group at 
the scale Mjj and hence does not affect the RGE. The (15,1,1) component 
of the 45-plet field H breaks Gps and so contributes to the RGE between the 
scale M c and My. The color singlet part of the (10,1,3) component of the 
126-plet field breaks the Glr symmetry. Then by extended survival hy- 
pothesis the color singlet component contribute to the RGE between the scale 
Mr to M c and all the components (10,1,3) contribute between the scales 



M c and My On the other hand our proposed mechanism, allows all of 
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the components of the field = (10,3,1) to remain light and contribute 
to the RGE at all energies between M\y and My. The bidoublet (1,2,2) 
field breaks the electroweak symmetry group and contributes at all energies 
to the RGE . For the correct quark-lepton mass relation we also require the 
bidoublet color singlet component of (15,2,2) to acquire a vev at the elec- 
troweak scale. For the potential we have, the (6,2,2) field will mix with the 
color triplet and anti-triplet components of the (15,2,2) field and one com- 
bination of these color triplet and anti-triplet fields will remain light. The 
color singlet and one combination of the triplet and anti-triplet component 
will then contribute to the RGE at all energies, while the other combination 
of the color triplet and anti-triplet will become heavy and contribute only 
between the energies M c and My. 



My -> M c 


M C ^M R 


Mr -> M W 


(1,1,1) 






(15,1,1) 






(10,1,3) 


(1,1,3,^1) 




(10,3,1) 


(6,3,l,- A /|) + (3,3,l, v /i) +(1,3,1,71) 


(6,3,-^/f) +(3,3, I v /|) + (l,3, v /|) 


(1,2,2) 


(1,2,2,0) 


(l,2,y§)+(l,2, =±y/l) 


(15,2,2) 


(3,2,2, v /|) + (3,2,2,- V / |) + (l,2,2,0) 


(3,2,| V / |) + (3,2,| v /|) + (l,2, yi) 
(3,2, ^V|)+(3,2, ^S + (l,2, 


(6,2,2) 







Table 1: Higgs scalars at various symmetry breaking scales. The U(l) quan- 
tum numbers are normalized from their embedding in S'O(IO). 

The normalization of the U(l) quantum numbers at the right handed 
breaking scale is fixed by the relation, 
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Combining results of Eq. [12] and Table. |T] it is easy to write down the 
following explicit form of the beta functions that regulate the evolutionary 
behaviour of the gauge couplings at various energy scales. We assume that 
the number of fermion generations is three. 



Mjj — 


>M C 


M c -> Mr 


Mr- 


>M Z 




11/3 


b 3cr = -29/6 


b3crw 


-29/6 


i>2Luc = 


11/3 


b 2 Lcr= 4/3 


b^hrw = 


= 4/3 


i>2Ruc= 


11/3 


b>2Rcr= 3/2 


blYrw = 


49/6 






blB-Lcr= 13 







Table 2: The modified beta functions {b N = inbjy) for the various groups 
at different energy scales. In the table we use the notation b^xy, where N 
represents the group (N for SU(N) and IS for U(l)s) and xy means the 
beta functions within the scales M x and M y . 

We consider both symmetry breaking scales, Mu and M c , to be very large 
so that Planck scale effects are not negligible. We start with the renormaliz- 
able 5*0 (10) invariant lagrangian, 



--Tr(F^F^) 



(15) 



and then include the non-renormalizable higher dimensional terms which 
have their origin in Planck scale physics. We consider only terms of dimension 
5 and 6, given by 

-Tr(F^F^) (16) 



2M. 



pi 
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1 1 

'2 Ml 



rj^ {Tr (F^F^) + Tr(i^0F^0)} + 



^ 2) Tr(0 2 )Tr(F^F^) + ?^Tr(F^0)Tr(F^ 



(17) 



where are dimensional couplings of the higher dimensional operators. 
When any Higgs scalar <p acquires vev <p , these operators induce effective 
dimension 4 terms modifying the boundary conditions at the scale O - 

The symmetry breaking at M v shifts the boundary condition of the SU{A) 
coupling constant with respect to the SU (2) couplings whereas the vevs of 
the 45-plet field H contribute to the relative couplings of the SU (3) and the 
U(l) constants. The Gps invariant effective lagrangian, modified by these 
higher dimensional operators, is given by, 

_ 1 (i + q) Tr (F£) F^) - 1(1 + e 2 ) Tr(if^) F^) 



;i + e 2 ) Tr(F(f) F^"") 



(18) 



where, 



e 4 = e (1) + e 



(2) _, i J2) 



and 



£2 = 



» 



.La) + L(2) 

2 4 a 



+ 2^ 



(2) 



[J 1 1 


2 M[/l 


. 1 257Tq;g - 





(n) 



Then the usual Gps lagrangian can be recovered with the modified coupling 
constants, 



gl(Mu) 

9 2 2r(Mu) 



^ 4 2 (M {7 )(l + e 4 )- 1 
^ 2 2 L (M [/ )(l + e 2 )- 1 
^(AfoXl + ea)- 1 



(19) 
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where the cji are the coupling constants in the absence of the nonrenormaliz- 
able terms and gi are the physical coupling constants that evolve below M\j. 
The modified boundary condition then reads, 

gl(Mu)(l + e 4 ) = gl L {Mu){l + e 2 ) = cfc R {M v ){\ + e 2 ) = g 2 . (20) 

At M c the symmetry group SU(4) C breaks down to SU(3) C x U(\)b-l 
when the (15,1,1) component of the 45-plet of Higgs field H acquires a vev. 
The SU(3) C x U(1)b-l invariant lagrangian is given by, 

_ 1 (i + 4) Tr(i^) F^") - 1(1 + ei) Tr(F« F^"") 



where, 



f ' - f '(2) , 19/(2) 

7e , i 2 )+12e' (2) + 126 , ; 2 \ 



and 



'(2) 



24A/|, 



207T« 4 



Mr 



! 2' 



M, 



/(2) 



where, i — a,b, c. Then the boundary condition at M c becomes, 



9 2 i(b-l)(M c )(1 



g 2 3c (M c )(l + e' 3 )=gl(M c ). 



The matching conditions at the scale Mr are not modified by the Planck 
scale effects and are given by, 



R) 



r) 



92r(Mr) + -g UB _ L) (M R ) 
r) 



5 

9m{M< 



(21) 



Using the above boundary conditions and the one loop renormalization 
group equation the unification coupling ajj can be related to the three cou- 
plings at the W mass scale {My/) through the following relations [15| (we 
have defined m,- 



ln^l 
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a~\M w ) = aa 1 (1 + 3/5e 2 + 2/5e 4 + 2/5(e'(l + e 4 ))) + (6/5 b 2ruc + 4/5(1 + e[) b 4uc ) m t 

+ (6/5 b 2rcr + 4/5 6 
a 2 ~ 1 (M w ) = otQ 1 (1 + e 2 ) + 2 b 2 i uc m uc + 2 6 2 Zcr rn cr + 2 & 2/riu m r 
a^ 1 (M w ) = a v l (1 + e 4 ) + 2 6 4uc m uc + 2 6 3cr m cr + 2 6 3cru) m r 



We define, 



and 



A = ay X [M w )-a- 1 [M z ) 



B = a 2 l{M w ) + ^a Y \M z ) - -a£(M w ) 

and relate them to the experimental numbers through the following equa- 
tions, 

sm e w = - - 



8 8 



1--— = aB. (22) 
8 a, V ; 

If we now take all the e s to be zero, then we have 

m„ = — 36.7 + 5.7m uc 
m cr = 147.3 — 11. 2m uc . 

There is no solution with positive m rw and m cr for any value of M uc with 
the constraints, sin 2 6w = 0.2334, a s = 0.12, and with the unification scale 
below the Planck scale . In other words this means that if we do not consider 
the effect of gravity, then in the presence of so many light Higgs scalars it is 
not possible to have unification of the gauge coupling constants. Thus if the 
three lepton decay mode of the proton is the explanation of the atmospheric 
neutrino problem, gravity effects modify the low energy predictions of the 
grand unified theories. 
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We now consider the effects of the Planck scale. For several choices of 
the e s it is possible to have a mass scale solution which may explain the 
atmospheric neutrino problem. To demonstrate this we present a few rep- 
resentative solutions in table 2. The light scalars contribute from the scale 
M z ~ 100 GeV in the RG equation. The unification scale Mu and the 
Gps breaking scale M c are very close to each other (i.e., m uc = hi^ is 
very small) and is considered here to be around Mjj ~ M c ~ 10 18 GeV 
(i.e., m rw + m cr + m uc = 39); the right handed breaking scale is around 
Mr ~ 10 13 GeV (i.e., m rw ~ 28). With these values of the mass scales the 
three lepton decay mode of the proton would be the most dominant decay 
mode (with t(P — > e + vv) ~ 10 31 yrs). Since the unification scale is quite 
high now, conventional proton decay modes are very much suppressed (with 
t(P — > e + 7r°) ~ 10 39 yrs.). Thus even though the three lepton decay mode 
of the proton explains the atmospheric neutrino anomaly there is no conflict 
with the non-observation of proton decay in other experiments. 



A 


62 


4 


e 4 




771 rw 


m cr 




-l 


-.75 


1 


-1 


55 


28.32 


9.73 


.93 


-l 


-.75 


1 


-1 


57 


28.04 


10.34 


.61 


-l 


-.75 


1 


-1 


59 


27.75 


10.95 


.28 


-.99 


-.75 


1 


-1 


59 


27.76 


10.95 


.28 


-.90 


-.75 


1 


-1 


59 


27.76 


10.94 


.28 


-.80 


-.75 


1 


-1 


59 


27.77 


10.93 


.28 



Table 3: Allowed ranges of parameters for unification 
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6 Conclusions 



We have presented an extension of the left-right symmetric model where the 
most dominant proton decay mode is through three leptons. The lifetime for 
this decay mode is large enough to explain the atmospheric neutrino anomaly. 
We have minimized the complete potential to check the consistency of the 
model. In the end we have carried out a renormalization group analysis 
to estimate the mass scales of the model. We have shown that when the 
gravity induced effects coming from the Plank scale physics are included in 
the renormalization group analysis, the couplings unify; an estimation of the 
several mass scales of the model then becomes possible. 
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Figure 1: Diagram giving P — > e^v^VL 
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